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Abstract — Convex relaxations of the optimal finger selection 
algorithm are proposed for a minimum mean square error 
(MMSE) Rake receiver in an impulse radio ultra-wideband 
system. First, the optimal linger selection problem is formu- 
lated as an integer programming problem with a non-convex 
objective function. Then, the objective function is approximated 
by a convex function and the integer programming problem 
is solved by means of constraint relaxation techniques. The 
proposed algorithms are suboptimal due to the approximate 
objective function and the constraint relaxation steps. However, 
they can be used in conjunction with the conventional finger 
selection algorithm, which is suboptimal on its own since it 
ignores the correlation between multipath components, to obtain 
performances reasonably close to that of the optimal scheme 
that cannot be implemented in practice due to its complexity. 
The proposed algorithms leverage convexity of the optimization 
problem formulations, which is the watershed between 'easy' and 
'difficult' optimization problems. 

Index Terms — Ultra-wideband (UWB), impulse radio (IR), 
MMSE Rake receiver, convex optimization, integer programming. 



I. Introduction 

Since the US Federal Communications Commission (FCC) 
approved the limited use of ultra-wideband (UWB) technology 
[1], communications systems that employ UWB signals have 
drawn considerable attention. A UWB signal is defined to be 
one that possesses an absolute bandwidth larger than 500MHz 
or a relative bandwidth larger than 20% and can coexist with 
incumbent systems in the same frequency range due to its 
large spreading factor and low power spectral density. UWB 
technology holds great promise for a variety of applications 
such as short-range high-speed data transmission and precise 
location estimation. 

Commonly, impulse radio (IR) systems, which transmit very 
short pulses with a low duty cycle, are employed to implement 
UWB systems ([2]-[6]). In an IR system, a train of pulses is 
sent and information is usually conveyed by the position or 
the polarity of the pulses, which correspond to Pulse Position 
Modulation (PPM) and Binary Phase Shift Keying (BPSK), 
respectively. In order to prevent catastrophic collisions among 
different users and thus provide robustness against multiple- 
access interference, each information symbol is represented by 
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a sequence of pulses; the positions of the pulses within that 
sequence are determined by a pseudo-random time-hopping 
(TH) sequence specific to each user [2]. The number Nf 
of pulses representing one information symbol can also be 
interpreted as pulse combining gain. 

Commonly, users in an IR-UWB system employ Rake re- 
ceivers to collect energy from different multipath components. 
A Rake receiver combining all the paths of the incoming 
signal is called an all-Rake (ARake) receiver Since a UWB 
signal has a very wide bandwidth, the number of resolvable 
multipath components is usually very large. Hence, an ARake 
receiver is not implemented in practice due to its complexity. 
However, it serves as a benchmark for the performance of 
more practical Rake receivers. A feasible implementation of 
multipath diversity combining can be obtained by a selective- 
Rake (SRake) receiver, which combines the M best, out of 
L, multipath components [7]. Those AI best components are 
determined by a finger selection algorithm. For a maximal 
ratio combining (MRC) Rake receiver, the paths with highest 
signal-to-noise ratios (SNRs) are selected, which is an optimal 
scheme in the absence of interfering users and inter-symbol 
interference (ISI). For a minimum mean square error (MMSE) 
Rake receiver, the "conventional" finger selection algorithm 
is to choose the paths with highest signal-to-interference- 
plus-noise ratios (SINRs). This conventional scheme is not 
necessarily optimal since it ignores the correlation of the 
noise terms at different multipath components. In other words, 
choosing the paths with highest SINRs does not necessarily 
maximizes the overall SINR of the system. 

In this paper, we formulate the optimal MMSE SRake as a 
nonconvex, integer-constrained optimization, where the aim is 
to choose the finger locations of the receiver so as to maximize 
the overall SINR. While computing the optimal finger selec- 
tion is NP-hard, we present several relaxation methods to turn 
the (approximate) problem into convex optimization problems 
that can be very efficiently solved by interior-point methods, 
which are polynomial time in the worst case, and very fast 
in practice. These optimal finger selection relaxations produce 
significantly higher average SINR than the conventional one 
that ignores the correlations, and represent a numerically 
efficient way to strike a balance between SINR optimality and 
computational tractability. 

The remainder of the paper is organized as follows. Sec- 
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Fig. 1. An example time-hopping impulse radio signal with pulse-based 
polarity randomization, where Nj = 6, Nc = 4, the time hopping sequence 
is {2,1,2,3,1,0} and the polarity codes are {-i-1,-h1,-1,-i-1,-1,-h1}. 
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tion II describes the transmitted and received signal models 
in a multiuser frequency-selective environment. The finger 
selection problem is formulated and the optimal algorithm 
is described in Section III, which is followed by a brief 
description of the conventional algorithm in Section IV. In 
Section V, two convex relaxations of the optimal finger se- 
lection algorithm, based on an approximate SINR expression 
and integer constraint relaxation techniques, are proposed. 
The simulation results are presented in Section VI, and the 
concluding remarks are made in the last section. 

II. Signal Model 

We consider a synchronous, binary phase shift keyed TH- 
IR system with K users, in which the transmitted signal from 
user k is represented by: 




(1) 

where ptxit) is the transmitted UWB pulse, Ek is the bit 
energy of user k, Tf is the "frame" time, Nf is the number of 
pulses representing one information symbol, and ^[j/jVj.j ^ 
{+1,-1} is the binary information symbol transmitted by 
user k. In order to allow the channel to be shared by many 
users and avoid catastrophic collisions, a time-hopping (TH) 



sequence {c^ '}, where ' G {0, 1, A^c — 1}, is assigned 
to each user. This TH sequence provides an additional time 
shift of Cj'^Tc seconds to the jth pulse of the fcth user where 
Tc is the chip interval and is chosen to satisfy Tc < Tf/Nc 
in order to prevent the pulses from overlapping. We assume 
Tf — NcTc without loss of generality. The random polarity 

(k) 

codes d)j are binary random variables taking values ±1 with 
equal probabiHty ([8]-[10]). 

Consider the discrete presentation of the channel, a^'^^ — 
[a^'^' • • • a^^^] for user fc, where L is assumed to be the number 
of multipath components for each user, and Tc is the multipath 
resolution. Then, the received signal can be expressed as 
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Fig. 2. The receiver structure. There are M multipath components, which 
are combined by the MMSE combiner. 



of the antenna, and n{t) is zero mean white Gaussian noise 
with unit spectral density. 

We assume that the time-hopping sequence is constrained 
to the set {0, 1, . . . , TVt - 1}, where Nt < Nc - L, so that 
there is no inter- frame interference (IFI). 

Due to the high resolution of UWB signals, chip-rate and 
frame rate sampling are not very practical for such systems. 
In order to have a lower sampling rate, the received signal can 
be correlated with template signals which enable symbol rate 
sampling of the output [11]. The template signal for the Ith 
path of the incoming signal can be expressed as 

^(t-jTf-cY>Tc-{l-l)Tc), 
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for the ith information symbol, where we consider user 1 
without loss of generality. In other words, by using a correlator 
for each multipath component that we want to combine, we 
can just use symbol rate sampling at each branch, as shown 
in Figure 12] 

Note that the use of such template signals results in equal 
gain combining (EGC) of different frame components. This 
may not be optimal under some conditions (see [12] for 
(sub)optimal schemes). However, it is very practical since it 
facilitates symbol-rate sampling. Since we consider a system 
that employs template signals of the form (|3}, i e. EGC of 
frame components, it is sufficient to consider the problem of 
selection of the optimal paths just for one frame. Hence, we 
assume Nf = 1 without loss of generality. 

Let £ = {li, . . . , denote the set of multipath compo- 
nents that the receiver collects (Figure |2j. At each branch, 
the signal is effectively passed through a matched filter (MF) 
matched to the related template signal in (|3} and sampled once 
for each symbol. Then, the discrete signal for the kh path can 
be expressed, for the ith information symbol, as^ 

ri = sfAhi + ni, 



(4) 



Tc-il- l)Tc) + ann{t), for = h,- -j hi, where A_= diag{V^ 



(2) [b[ 



and 71/ ~ A/'(0 , a'^). s/ is a X x 1 vector. 



where p^x (t) is the received unit-energy UWB pulse, which is 
usually modelled as the derivative of ptx (t) due to the effects 



-Note that the dependence of r; on the index of the information symbol, 
i, is not shown explicitly. 



which can be expressed as a sum of the desired signal part 
(SP) and muhiple-access interference (MAI) terms: 



(SP) 
Si = s\ 
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where the fcth elements can be expressed as 
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and 



fc= 1 



(5) 
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with /j „j being the indicator function that is equal to 1 if the 
TOth path of user k collides with the Ith path of user 1, and 
otherwise. 



III. Problem Formulation and Optimal Solution 



The problem is to choose the optimal set of multipath 
components, C = {li, . . . , Im}, that minimizes the bit error 
probability (BEP) of the system. In other words, we need 
to choose the best samples from the L received samples r/, 
I = 1, . . . , L, as shown in @. 

To reformulate this combinatorial problem, we first define 
an M X L selection matrix X as follows: M of the columns 
of X are the unit vectors ei, . . . ,eM (e^ having a 1 at its 
ith position and zero elements for all other entries), and the 
other columns are all zero vectors. The column indices of the 
unit vectors determine the subset of the multipath components 
that are selected. For example, for L — A and M — 2, X- — 

Q Q ^ Q chooses the second and third multipath components. 

Using the selection matrix X, we can express the vector of 
received samples from M multipath components as 

r = XSAb,; + Xn, (8) 

where n is the vector of thermal noise components n = 
[ni ■ ■ ■ ni,]'^ , and S is the signature matrix given by S = 
[si • • -sl]"^, with s; as in (|5}- 
Using (|5}-(Q> ll} can be expressed as 

r = feW V^Xa^i) + XS(*^^^) Ab, + Xn, (9) 

where S^*^"^^) is the MAI part of the signature matrix S. 
Then, the linear MMSE receiver can be expressed as 

h = sign{0^r}, (10) 

where the MMSE weight vector is given by [13] 



(11) 



with R being the correlation matrix of the noise term: 



The SINR of the system can be expressed as 
SINR{X) = :^(a(i))^X^ 

I+^XS(*^-4^)A2(S(^^^)fX^) XaW. (13) 



Hence, the optimal path/finger selection problem can be for- 
mulated as 



maximize SINR(X.), 



(14) 



where X has the previously defined structure. 

Note that the objective function to be maximized is not 
concave and the optimization variable X takes binary values, 
with the previously defined structure. In other words, two 
major difficulties arise in solving il4i globally: nonconvex 
optimization and integer constraints. Either makes the problem 
NP-hard. Therefore, it is an intractable optimization problem 
in this general form. 

IV. Conventional Algorithm 

Instead of the solving the problem in il4\ . the "conven- 
tional" finger selection algorithm chooses the M paths with 
largest individual SINRs, where the SINR for the Ith path can 
be expressed as 



SINRi 



(MAI).^j,j^2^(MAI) 



(15) 



R = XS(^'^^^)a2(S(^'^^^))^X^ + all. (12) 



forZ = 

This algorithm is not optimal because it ignores the corre- 
lation of the noise components of different paths. Therefore, 
it does not always maximize the overall SINR of the system 
given in ( I13> . For example, the contribution of two highly 
correlated strong paths to the overall SINR might be worse 
than the contribution of one strong and one relatively weaker, 
but uncorrelated, paths. The correlation between the multipath 
components is the result of the MAI from the other users in 
the system. 

V. Relaxations of Optimal Finger Selection 

Since the optimal solution in il4l is quite difficult, we first 
consider an approximation of the objective function in il3i . 
When the eigenvalues of ^XS(*^'4^) A2(S(*^'4^))^X^ are 
considerably smaller than 1, which occurs when the MAI 
is not very strong compared to the thermal noise, we can 
approximate the SINR expression in ( I13t as follows"^: 

\ - ^XS(*^^^)A2(S(*^^^))^X^j Xa(i), (16) 

^More accurate approximations can be o btain ed by using higher order 
series expansions for the matrix inverse in 1131 . However, the solution of 
the optimization problem does not lend itself to low complexity solutions in 
those cases. 



which can be expressed as 
SINR{X) « ^{ (a(i))^X^Xa(i) 

- ±a(i)x^XS(*^^^)A2(S(*^^^))^X^Xa(i) }. (17) 

Note that the approximate SINR expression depends on X 
only through X-^X. Defining x = [xi • • • ccl]^ as the diagonal 
elements of X^X, x = diagjX^X}, we have = 1 if the 
ith path is selected, and Xi ~ otherwise; and J2i=i — 
Then, we obtain, after some manipulation. 



SINRM = ^ 



T ^ 

q X 



x^ Px 



(18) 



where q = [{al^f---{a'^'YY' and P 
diag{ai^) • • •ai^^}S(^^^^)A2(S(*^^^))^diag{a« • ••««}. 

Then, we can formulate the finger selection problem as 
follows: 



mmimize 



-i-x"^Px — x^'^q 



subject to x"'"! = M, 
Xi e {0,1}, 



1, 



(19) 



Note that the objective function is convex since P is positive 
definite, and that the first constraint is linear. However, the 
integer constraint increases the complexity of the problem. 
The common way to approximate the solution of an integer 
constraint problem is to use constraint relaxation. Then, the 
optimizer will be a continuous value instead of being binary 
and the problem ( I19> will be convex. Over the past decade, 
both powerful theory and efficient numerical algorithms have 
been developed for nonlinear convex optimization. It is now 
recognized that the watershed between "easy" and "difficult" 
optimization problems is not linearity but convexity. For 
example, the interior-point algorithms for nonlinear convex 
optimization are highly efficient, both in worst case complexity 
(provably polynomial time) and in practice (very fast even for 
a large number variables and constraints) [14]. Interior-point 
methods solve convex optimization problems with inequality 
constraints by applying Newton's method to a sequence of 
equality constrained problems, where the Newton's method is 
a kind of descent algorithm with the descent direction given 
by the Newton step [14]. 

We consider two different relaxation techniques in the 
following subsections. 

A. Case-1: Relaxation to Sphere 

Consider the relaxation of the integer constraint in il9\ to 
a sphere that passes through all possible integer values. Then, 
the relaxed problem becomes 

minimize ^x"^Px — x^q 



subject to X 1 = M, 

(2x- l)^(2x- 1) < L. 



Note that the problem becomes a convex quadratically con- 
strained quadratic programming (QCQP) [14]. Hence it can 
be solved for global optimality using interior-point algorithms 
in polynomial time. 

B. Case-2: Relaxation to Hypercube 

As an alternative approach, we can relax the integer con- 
straint in ( I19> to a hypercube constraint and get 

1 

mmirmze — 

subject to x^l — M, 

xe [0,1]^, 



-x^Px — x^q 



(21) 

where the hypercube constraint can be expressed as x )^ 
and X ^ 1, with y ^ z meaning that yi > zi, . . . ,yL > zl- 
Note that the problem is now a linearly constrained quadratic 
programming (LCQP), and can be solved by interior-point 
algorithms [14] for the optimizer x*. 

C. Dual Methods 

We can also consider the dual problems. For the relaxation 
to the sphere considered in Section IV-AI the Lagrangian for 
(I20> can be obtained as 

/:(x, A, v) = x^ ^^P + 4iylj X - x^(q - Al + Avl) - MX, 

(22) 

where \ E TZ and ly E TZ^. 

After some manipulation, the Lagrange dual function can 
be expressed as 

-1 



i[q+(A + 4i.)lFf^P + 4i.I 



[q+ (A + 4i/)l] - A/A, 
Then, the dual problem becomes 



(23) 



mmimize 

i [q + (A + Aiy)lf (^^P + Aiylj [q + (A + 4t/)l] + MX 

(24) 
(25) 



subject to > 0, 



(20) 



which can be solved for optimal A and ly by interior point 
methods. Or, more simply, the unconstraint problem i24l can 
be solved using gradient descent algorithm, and then the 
optimizer D is mapped to v* — max{0, D}. 

After solving for optimal A and /i, the optimizer x* is 
obtained as 

X* = i (^^P + 4i.*lj [q+(A*+4;.*)l]. (26) 

Note that the dual problem ( I24> has two variables, A and v, 
to optimize, compared to L variables, the components of x, 
in the primal problem i2Q\ . However, an L x L matrix needs 
to be inverted for each iteration of the optimization of ( I24t . 
Therefore, the primal problem can be preferred over the dual 
problem in this case. 



Similarly, the dual problem for the relaxation in Section lV-BI 
can be obtained from i21\ as 



minimize 

^2 



-{(1 + H u Al)^p-^(q + /X - - Al) + MA 



subject to fj,,v y 0. 



(27) 
(28) 



It is observed from ( I27> that there are 2L + 1 variables and 
also L X L matrix inversion operations for the solution of the 
dual problem. Therefore, the simpler primal problem J21> is 
considered in the simulations. 

D. Selection of Finger Locations 

After solving the approximate problem il9i by means 
of integer relaxation techniques mentioned above, the finger 
location estimations are obtained by the indices of the M 
largest elements of the optimizer x*. 

Both the approximation of the SINR expression by J16t and 
the integer relaxation steps result in the suboptimality of the 
solution. Therefore, it may not be very close to the optimal 
solution in some cases. However, it is expected to perform 
better than the conventional algorithm most of the time, since 
it considers the correlation between the multipath components. 
However, it is not guaranteed that the algorithms based on 
the convex relaxations of optimal finger selection always beat 
the conventional one. Since the conventional algorithm is 
very easy to implement, we can consider a hybrid algorithm 
where the final estimate of the convex relaxation algorithm 
is compared with that of the conventional one and the one 
that minimizes the exact SINR expression in ( I13t is chosen as 
the final estimate. In this way, the resulting hybrid suboptimal 
algorithm can get closer to the optimal solution. 

VI. Simulation Results 

The simulation results are performed to evaluate the per- 
formance of different finger selection algorithms for an IR- 
UWB system with = 20 and Nf ^ 1. There are 5 
equal energy users in the system (K — 5) and the users' 
TH and polarity codes are randomly generated. We model the 
channel coefficients as ai = sign(a;)|Q!;| for I = 1,...,L, 
where sign (a;) is ±1 with equal probability and |a;| is 
distributed lognormally as £Af{tii,a^). Also the energy of 
the taps is exponentially decaying as E{|aip} = D,oe~^^^^^\ 
where A is the decay factor and X^^^i ^{1*^'^} ~ 1 (^^ 



ilo = (1 



^)/(l 



-XL 



we have A — 0.1, a 
t^i = 0.5 iH^Ef^ 



)). For the channel parameters, 
= 0.5 and fii can be calculated from 
- A(; - 1) - 2cr2j , for / = 1,...,L. 
We average the overall SINR of the system over different 
realizations of channel coefficients, TH and polarity codes of 
the users. 

In Figure |3] we plot the average SINR of the system for 
different noise variances when M — 5 fingers are to be chosen 
out of L = 15 multipath components. As is observed from 
the figure, the convex relaxations of optimal finger selection 
result in SINR values reasonable close to those of the optimal 




Fig. 3. Average SINR versus Eh /No for A4 — 5 fingers. The channel 
has L — 15 multipath components and the taps are exponentially 
decaying. The IR-UWB system has Nc = 20 chips per frame and 
Nf — 1 frame per symbol. There are 5 equal energy users in the 
system and random TH and polarity codes are used. 



exhaustive search scheme. Note that the gain by using the 
proposed algorithms over the conventional one increases as 
the thermal noise decreases. This is because when the thermal 
noise get less significant, the MAI becomes dominant, and 
the conventional technique gets worse since it ignores the 
correlation between the MAI noise terms when choosing the 
fingers. 

Next, we plot SINR of the proposed suboptimal and con- 
ventional techniques for different finger numbers in Figure 
0] where there are 50 multipath components and Ei,/Nq — 
20. The number of chips per frame, Nc, is set to 75, and 
all other parameters are kept the same. In this case, the 
optimal algorithm takes a very long time to simulate since it 
needs to perform exhaustive search over many different finger 
combinations (therefore not implemented). The improvement 
using convex relaxations of optimal finger selection over 
the conventional technique decreases as M gets large since 
the channel is exponentially decaying and the most of the 
significant multipath components are already combined by all 
the algorithms. 

Finally, we consider a MAI-limited scenario, where there 
are 10 users with Ei = 1 and Ek = 10 Vfc, and all the 
parameters are as in the previous case. Then, as shown in 
Figure |5] the improvement by using the suboptimal finger 
selection algorithms increase significantly. The main reason 
for this is that the suboptimal algorithms consider, although 
approximately, the correlation caused by MAI whereas the 
conventional scheme simply ignores that. 

VII. Concluding Remarks 

Optimal and suboptimal finger selection algorithms for 
MMSE-SRake receivers in an IR-UWB system are consid- 
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Fig. 4. Average SINR versus number of fingers M, for Et/No — 
20dB, A^c = 75 and L — 50. All the other parameters are the same 
as those for Figure |3| 
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Fig. 5. Average SINR versus number of fingers AI. There are 10 
users with each interferer having lOdB more power than the desired 
user. All the other parameters are the same as those for Figure |4] 



ered. Since UWB systems have large number of multipath 
components, only a subset of those components can be used 
due to complexity constraints. Therefore, the selection of the 
optimal subset of multipath components is important for the 
performance of the receiver. We have shown that the optimal 
solution to this finger selection problem requires exhaustive 
search which becomes prohibitive for UWB systems. There- 
fore, we have proposed approximate solutions of the problem 
based on the Taylor series approximation and integer constraint 
relaxations. Using two different integer relaxation approaches, 
we have introduced two convex relaxations of the optimal 



finger selection algorithm. Implementing these suboptimal 
algorithms on top of the conventional scheme, we can get 
close to the optimal solution, with much lower complexity. 

The two contributions of the paper are the formulation of the 
optimal problem and the convex relaxations. In the first, the 
formulation is globally optimal but the solution methods for 
non-convex nonlinear integer constrained optimization must 
use heuristics to get to a locally optimal solution because oth- 
erwise computational load for global optimality is too much. 
In the second, the formulation is relaxed, but the interior-point 
methods efficiently computes the globally optimal solution for 
these relaxations. 
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